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Abstract. - The upper Vienne River is the westernmost drainage in the natural distribution of European grayling 
Thymallus thymallus (Linnaeus, 1758). The grayling extension across this sub-basin of the Loire catchment has 
been dramatically reduced through habitat degradation, which began in the 19 th century with the construction 
of milldams and eventually large hydroelectric dams, which flooded former grayling habitat (i.e. fast flowing 
waters) under increasingly larger reservoirs. Throughout the 20 th century, habitat degradation, fragmentation and 
pollution have strongly affected the remaining isolated populations. Since the 1960s, fishing authorities have 
attempted to compensate for these declines with stocking programs. Most stocked fish stem from imported eggs 
or fry of unknown origin. However, the efficacy of this policy has never been assessed. Here, we aim to assess 
the genetic traces of these stocking programs as an indirect measure of the long-term survival of stocked fish. 
Sampled grayling from four localities in the Upper Vienne basin were analysed at both mtDNA (N = 37, control 
region) and nDNA levels (N = 88,10 microsatellites), and compared to populations representative of surrounding 
drainage basins or fish farms. Microsatellite analysis demonstrated low genetic diversity within samples, typical 
for the Loire drainage. Both Factorial Correspondence and “Structure” analyses presented Upper Vienne gray¬ 
ling as a comprehensive genetic unit distinct from diverse samples from the Rhine, Rhone and Danube basins, as 
well as other regions of the Loire catchment and the sampled hatcheries. Within the Vienne, only the uppermost 
sample differed from the other sites. Among the 37 individuals sequenced from the Upper Vienne, only three 
control region haplotypes were identified, all closely related to other haplotypes observed in Loire and Allier (a 
Loire tributary) rivers. These results strongly support that grayling in the Upper Vienne drainage represent native 
genetic stocks with no signs of non-native introgression. The only evidence potentially supporting a stocking 
contribution is the co-occurrence of a haplotype in the Allier catchment and Vienne River, although we cannot 
exclude that this haplotype is also native to both sites. These results suggest that five decades of stocking have 
provided no added reproductive value to the Vienne populations. We thus hypothesize that long-term survival of 
stocked fish has either been low or non-existent. 


Resume. - Integrite genetique des populations d’ombre commun ( Thymallus thymallus ) du bassin de la Vienne 
apres cinq decennies d’alevinages. 

Le bassin amont de la Vienne constitue la limite occidentale de l’aire de repartition naturelle de Tombre 
commun Thymallus thymallus (Linnaeus, 1758) en Europe, mais depuis le xix e siecle la construction de multi¬ 
ples moulins puis de grands barrages hydroelectriques a singulierement reduit et fragments le lineaire de cours 
d’eau habitable par cette espece rheophile. De plus, les populations subsistantes ont ete exposees aux pollutions 
industrielles des xix e et xx e siecles, ce qui a conduit a Textinction de beaucoup d’entre elles. Pour tenter de com- 
penser ce declin, des le milieu des annees soixante les autorites locales de la peche ont procede regulierement a 
des deversements d’alevins. Ces derniers ont ete fournis par divers pisciculteurs publics ou prives important leurs 
oeufs de l’exterieur, mais sans fournir dedications precises quant a leur origine exacte. Par ailleurs, Tefficacite 
de ces alevinages n’a jamais ete reellement evaluee. Le but de cette etude etait done d’examiner Timpact gene¬ 
tique de ces repeuplements. Quatre populations du bassin de la Vienne ont ete analysees a la fois au niveau de 
l’ADN nucleaire (10 micro satellites) et de TADN mitochondrial (region de controle) et comparees a des popu¬ 
lations representatives des bassins limitrophes ou des productions de piscicultures. Le sequengage de TADN 
mitochondrial de 37 individus “sauvages” n’a fourni que trois haplotypes, ce qui est peu mais comparable a ce 
qui est observe dans le bassin de la Loire et aussi souvent ailleurs. Ces trois haplotypes sont clairement apparen- 
tes aux haplotypes observes dans les bassins de TAllier et de la Loire en Auvergne, et par consequent eloignes 
de tous les autres haplotypes connus jusqu’a present en Europe. Par ailleurs, deux de ces haplotypes semblent 
specifiques du bassin de la Vienne tandis que le troisieme, observe seulement dans la population la plus aval de la 
Vienne, etait deja connu du haut Allier. L’analyse des microsatellites de 88 individus “sauvages” a egalement fait 
apparaitre une faible diversite dans toutes les populations, mais comparable a celle des autres populations lige- 
riennes. Que ce soit par Analyse Factorielle des Correspondances ou par Analyse “Structure”, les populations de 
la haute Vienne forment un ensemble homogene bien distinct de toutes les autres populations considerees. Tout 
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au plus la population amont de la Vienne sort-elle un peu du lot. Ces resultats prouvent que les populations de 
la haute Vienne sont bien natives de ce bassin, meme si elles s’apparentent aux autres populations naturelles du 
bassin de la Loire. Le seul indice d’une possible trace genetique d’alevinages anterieurs reside dans la presence 
chez quelques poissons de la Vienne aval d’un haplotype commun au haut Allier, mais on ne peut exclure que sa 
presence en Vienne y soit naturelle. L’absence d’autres traces de contributions genetiques allochtones a la Vienne 
suggere que cinq decennies d’alevinages en Ombres n’ont en rien contribue au maintien des populations locales. 
La survie ou le succes reproducteur des poissons de verses semblent avoir ete nuls, quels que soient la souche ou 
le stade de developpement employes. La presence d’une population autochtone semble un facteur determinant 
dans cet echec puisque de nombreux exemples d’implantations reussies sont a signaler dans des cours d’eau ou 
l’espece etait absente auparavant. 


For well over half a century, the viability and efficacy of 
stocking hatchery-reared salmonids into freshwaters have 
been a topic of both extensive research and controversy 
(Miller 1953, 1958; Vincent, 1987; Weiss and Schmutz, 
1999 and references therein). Despite several reviews on the 
topic (e.g. Cresswell, 1981; Finnegan and Stevens, 2008; 
Hansen et al., 2009), generalizations on the overall success 
of stocking (from a management perspective) or the extent 
of its negative consequences have been elusive. For brown 
trout for instance, Naslund (1998) reported experimental 
survival rates between 0 and 26% for 0+, 1+ and 2+ old 
stocked fish after 12 months across several Swedish streams, 
and Weiss and Schmutz (1999) reported similarly disparate 


rates (<1% to 19%) for 2+ and 3+ old hatchery-reared fish in 
two Austrian streams, compared to much higher survival for 
their wild counterparts. Champigneulle et al. (2002) pointed 
out that “Atlantic” brown trout introduced as fingerlings in 
the Doubs River suffered a higher angling pressure than the 
native “Mediterranean” trout during the following years. As 
a result, larger trout found in the population were mostly 
“Mediterranean”. This trend might contribute to the genetic 
segregation observed in that river by Largiader and Scholl 
(1996). Harbicht et al. (2014), in a meta-analysis of 37 com¬ 
parisons of survival between wild, hybridized and hatchery- 
reared salmonids, concluded that wild fish almost always 
exhibit higher survival rates than hybridized or hatchery - 


Figure 1. - Distribution of grayling in 
Europe (green) and its former distri¬ 
bution in France around 1900 (red). 
Sample sites (white stars = wild popu¬ 
lations, black stars = hatcheries) are 
numbered following table II. 
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reared fish. Additional to estimating survival of hatchery- 
reared fish, the degree of introgression has been investigated 
across many different salmonid fishes. For instance, some 
authors suggest that the intensity of stocking correlates with 
non-native introgression (Marie et al., 2010), while oth¬ 
ers invoke environmental integrity as a potential barrier for 
introgression (Seehausen et al., 2008; Bittner et al., 2010; 
Winkler et al., 2011), as well as various biological isolating 
mechanisms (Winkler et al., 2011). In any case, the ration¬ 
ality of stocking programs as a means of compensating (or 
mitigating) pervasive habitat changes primarily caused by 
human activity remains questionable, especially in flowing 
waters. 

Studies on the success or effects of stocking European 
grayling Thymallus thymallus (L., 1758) are much less 
numerous as those for salmon and trout (but see Carlstein, 
1997; Turek et al., 2010), but presumably the efficacy of 
stocking for this species must be at least as problematic as 
it is for other salmonid species. European grayling is consid¬ 
ered a sensitive salmonid confined to the large, clean, cold, 
and fast running rivers of North and Central Europe (Fig. 1; 
Persat et al., 1978). Its requirements are so specific that Huet 
(1949) created a specific zone for it in his river zonation 
scheme. 

Nonetheless, and in spite of its increasing popularity 
as a sport fish, grayling populations have suffered serious 
declines, especially in Central Europe and along the southern 
limits of its distribution (see Weiss et al. (2013) and refer¬ 
ence therein; Persat (1996) for France). The causes of these 
declines are considered to be multi-factorial and include 
various river engineering activities, increased predation 
pressure by cormorants Phalacrocorax carbo (L., 1758), 



Figure 2. - Close-up of grayling distribution and four sampled sites 
(stars) within the upper Vienne district in 2012 (relief background 
from IGN-Geoportail). 


interspecific competition with introduced salmonids such 
as rainbow trout Oncorhynchus my kiss (Walbaum, 1792) 
and climate change (Uiblein et al., 2001). Due to the small 
size of their larvae, rearing grayling is generally considered 
more difficult than trout, charr and salmon, and thus grayling 
hatcheries in France, at least initially, rarely relied on rear¬ 
ing local strains, preferring imported eggs or fry from longer 
established sources in Central and Northern Europe, usually 
centred near large natural spawning grounds. Thus, com¬ 
mercially available grayling strains have been quite diverse 
across time and space compared to, for example, brown and 
rainbow trout, which often stem from centralized industrial 
production facilities. Furthermore, extensive genetic investi¬ 
gations on European grayling over the last 15-20 years clear¬ 
ly demonstrated that their natural genetic diversity, espe¬ 
cially at the lineage or basin level, far exceeds that of other 
European salmonids (Susnik et al ., 1999; Gross et al., 2001; 
Weiss et al., 2002, 2013; Gum et al., 2005,2009). Some of 
these investigations have clearly highlighted cross-basin 
introgression due to stocking activities, and a recent report 
covering regions of Northern Italy and Western Austria even 
has demonstrated a nearly complete local displacement of 
native stocks with hatchery strains (Meraner et al., 2014; 
Weiss et al., 2015). However, this same report highlighted 
nearly pure or pure genetic stocks in some systems where 
frequent management activities have taken place. This sup¬ 
ports, for example, the observations of Meraner et al. (2013) 
in the Upper Drau River in Austria, where a purely native 
sample of grayling was reported, despite the fact that the 
whole river has been extensively managed with mixed and 
foreign stocks of fish. Thus, there is no clear consensus on 
the long-term survival of stocked grayling, or on the levels 
of genetic introgression observed in grayling populations 
that have been stocked with foreign strains. 

In France, stocking of grayling started slowly in 1960s 
with eggs imported mainly from Central Europe, a practice 
that became more common after 1970. A progressive shift to 
local production of eggs from fish matured in captivity took 
place after 1980, and was mainly driven by fishing authori¬ 
ties (Conseil Superieur de la Peche, and Angling Federa¬ 
tions; Carmie et al., 1985). However, the fickle nature of the 
production limited these activities, and private fish farms 
using imported eggs still provide a significant amount of the 
fry stocked into French rivers. A first genetic-based evalu¬ 
ation of the potential contribution of stocking within the 
partitioned river network of the French Upper-Rhone River 
was attempted by Persat and Eppe (1997), using enzymatic 
markers. However, this effort primarily showed a significant 
difference between the main channel and tributaries, but this 
difference was not due to the local hatchery: indeed, almost 
no individuals from these regular stockings were found. 

Ongoing investigations on the genetic diversity of the 
native grayling populations in France using mtDNA and 
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microsatellite markers confirmed significant divergence 
between the three major drainages (Rhine, Rhone and Loire) 
and at a finer level, within each drainage according to their 
size and tributary complexity (Cattaneo et al., 2011; Persat 
and Weiss, unpubl. data). Thus, the natural genetic struc¬ 
ture of grayling in France makes it rather easy to distinguish 
native versus non-native genotypic profiles. This is espe¬ 
cially true for the Loire drainage (Fig. 1), as it is naturally 
inhabited by one of the most genetically divergent lineages 
in Europe (Weiss et al., 2002). 


Table I. - Amount of grayling stocked per year in the Upper Vienne 
drainage basin (mostly with fry before 1990, mainly with finger- 
lings after 1990). 


Year 

Amount 

Year 

Amount 

Year 

Amount 

1966 

10,000 

1982 

98,000 

1998 

? 

1967 

30,000 

1983 

85,000 

1999 

ca. 6,000 

1968 

40,000 

1984 

76,000 

2000 

ca. 4,500 

1969 

20,000 

1985 

74,000 

2001 

ca. 26,000 

1970 

60,000 

1986 

55,200 

2002 

16,000 

1971 

65,000 

1987 

43,000 

2003 

9,000 

1972 

120,000 

1988 

94,000 

2004 

20,000 

1973 

64,000 

1989 

? 

2005 

22,000 

1974 

100,000 

1990 

? 

2006 

14,000 

1975 

31,200 

1991 

? 

2007 

17,000 

1976 

40,000 

1992 

? 

2008 

6,300 

1977 

77,000 

1993 

? 

2009 

12,000 

1978 

55,000 

1994 

12,400 

2010 

20,000 

1979 

67,000 

1995 

7,000 

2011 

10,000 

1980 

110,920 

1996 

ca. 2000 

2012 

10,000 

1981 

93,500 

1997 

? 




Table II. - Samples used in this study, including the name of the population, the major drainage, 
geographic coordinates and the number of individuals typed for either micro satellites or the 
complete mtDNA control region. 


Population 

Drainage 

Pop. Code 

Coordinates 

N msats 

N mtDNA 

Vienne 






Combade 

Loire 

1 

1°37 , 40”E -45°42’47”N 

37 

11 

Vienne upstream 

Loire 

2 

1°37 , 22”E-45°45 , 45”N 

34 

11 

Vienne downstream 

Loire 

3 

1°25 , 08”E-45°52 , 55”N 

15 

13 

Maulde 

Loire 

4 

1°46 , 22”E-45°50 , 35 ,, N 

2 

2 

Allier (Vabres) 

Loire 

5 

2°4W26”E -44°55’05”N 

41 

10 

Upper Loire 

Loire 

6 

3°55 ’30”E-44°51 ’54”N 

43 

19 

Ance du Nord 

Loire 

7 

4°01 ’52”E-45° 15 ’44”N 

- 

7 

Ain (Mollon) 

Rhone 

8 

5°14 , 52”E-45°57 , 20”N 

42 

17 

Besse 

Hatchery 

9 

2°55 , 25”E-45°30 , 55”N 

29 

11 

Obenheim 

Hatchery 

10 

7°41 ’ 16”E-48°21 ’20”N 

40 

11 

Chauvey (Appenans) 

Hatchery 

11 

6°33 , 58”E-47°26 , 48 ,, N 

33 

7 

Chazey-Bons 

Hatchery 

12 

smow’e^^s^t’n 

43 

10 

Kohlenbach 

Danube 


12°23 , 00”E-47°39 , 56”N 

25 

15 

Total 

384 

144 


In order to test the genetic impact of stocking on Thymal- 
lus thymallus , we selected populations of the remote Upper 
Vienne River sub-basin, which is isolated from the other 
native populations of the Loire-Allier sub-basin by approx¬ 
imately 800 km of riverine distance. Moreover, the Upper 
Vienne is under jurisdictional control of a single angling 
federation, which was able to provide us an overview of the 
stocking operations performed over the last 50 years (Tab. I). 
The fish were provided by more than five different public or 
private fish farms across time, four of which we have sam¬ 
pled within this study. These sources are not exhaustive, as 
local angler societies may have contributed stocked fish on 
their own initiative. Mostly, the origin of the fish was unre¬ 
corded and differed across time and facilities. We gathered 
samples from four farms that still currently producing gray¬ 
ling to represent gene pools of potential relevance to histori¬ 
cal stocking, whereby only one (Besse farm) has been offi¬ 
cially involved in recent stockings into the Vienne basin. 


MATERIAL AND METHODS 

Samples from a total of 384 T. thymallus individuals 
from 13 locations, with a focus on the Loire basin, were 
included in the study (Figs 1,2; Tab. II). Samples included 
seven locations in the Loire catchment (four of which come 
from the Upper Vienne), four fish farms, and two refer¬ 
ence populations representing the Rhone drainage and the 
northern alpine Danube lineage (see Weiss et al., 2002 for 
an overview of major genetic lineages). Several individuals 
from the Ance du Nord River in the Loire basin were addi¬ 
tionally included (for mtDNA analysis only) as reference 
material. We also included a Finn¬ 
ish sequence because Finnish eggs 
were sometimes imported by local 
farms. 

For the Vienne basin, sampling 
was done by electric fishing in Nov. 
2011. The highest apparent den¬ 
sity (individuals collected per hour) 
was found in the Combade River, 
whereas the lowest density (two 
individuals collected over several 
hundred meters) was recorded 
in the Maulde River in the only 
short river reach still inhabited by 
grayling between two large dams 
(Fig. 2). Fish were aged based on 
scale analysis (most fish were 1+ 
or 2+ old, i.e. from 2010 or 2009 
generations) and each fish > 20 cm 
was photographed. Small fin-clips 
were taken from the tip of the pel- 
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vie fin and stored in 96% ethanol before releasing the fish 
unharmed into the river. 

Fish farm samples were generously provided by angling 
societies, and all represent potential source material for 
past or present stocking programs in the Upper Vienne 
basin. These are not necessarily representative of all facili¬ 
ties that have produced grayling from imported eggs, as 
many facilities no longer exist. The samples also represent 
“batches” and may not fully cover the genetic diversity of 
each facility - grayling egg batches in such farms are rou¬ 
tinely separated from one another, to avoid mixing healthy 
and unhealthy broods (Carmie et al., 1985) and thus differ¬ 
ent tanks may represent different gene pools. Three of the 
farms claim to produce their brood from local stocks; i.e. the 
Rhine River (Obenheim fish farm, Bas-Rhin Federation), the 
Allier River (Besse-en-Chandesse, Puy-de-Dome Federa¬ 
tion) and the Rhone River (Chazey-Bons, Ain Federation). 
The fourth farm (Appenans, Chauvey) is a private facility 
mostly dependent on imported eggs, though some produc¬ 
tion is claimed to stem from local Doubs (Rhone drainage) 
stocks. 


Genomic DNA was isolated from preserved tissue using 
a high salt (ammonium acetate) extraction protocol, modified 
from Miller et al. (1988). A total of 12 published microsatel¬ 
lite loci were combined into two multiplex reactions as listed 
in table III and applied to the screening of all 384 samples. 
PCR conditions were 5-30 ng template DNA, 2.5 pi 2x Type- 
it Multiplex Master Mix (Qiagen), 0.2 pM of each primer in 
a final reaction volume of 5 pi. The PCR cycles consisted of 
an initial denaturation step of 95°C for 5 min followed by 25 
cycles of denaturation at 95°C for 30 sec, annealing at 59°C 
(7-plex) or 60°C (5-plex) for 1 min 30 sec and extension at 
72°C for 30 sec and a final extension step at 60°C for 30 min. 
The alleles were scored using an internal size standard (Gen- 
eScan-500 ROX, Applied Biosystems) and analyzed with 
GeneMapper® software version 3.7 (ABI). 

Observed and expected heterozygosity, the mean number 
of alleles and deviations from Hardy-Weinberg Equilibrium 
(HWE) (using the intrapopulation fixation index Fis) were 
calculated using GENETIX version 4.05 (Belkhir et al., 
1996). Statistical significance was obtained with 10 000 per¬ 
mutations and a table-wide Bonferroni correction. To display 
genetic relationships among individuals we carried out a Fac- 


Table III. - Microsatellite locus description including the arrangement of the 7-plex and 5-plex reactions, the locus name and repeat motif, 
the primer sequence, the fluorescent-labelling dye, annealing temperature, original reference and focal species used to clone the locus. 
Note that the published names of Thy 1 (BFRO004), Thy62 (BFRO010), and Thy54 (BFRO018) differ from those given here. 



Locus 

Repeat 

motif 

Primer sequence (5"-3") 

Dye 

AT 

Reference 

Focal 

species 


Tth447 

(TG)n 

F: CTTGATTGCCATTGGATTGT 

R: GTTT CTT C A AC ATCCTT GTCGCCT CTA 

HEX 

59°C 

Junge et al. (2010) 

T. thymallus 


Tth305 

(CTTT) 8 

F: CTTT G A ATAT G ATGCGT G A AC 

R: GTTT CTT G AGTATACTGC AG ATAG ACC A 

NED 

59°C 

Junge et al. (2010) 

T. thymallus 

7 - Plex 

Thyl 

(CA)n 

F: CGCATCTGTATGAAAAACCT 

R: TGGTTTGGTAGGAGTTTCGT 

FAM 

59°C 

Susnik et al. (1999a) 

T. thymallus 

Tar 101 

(CTTT) 22 

F: CAGAGCACACCAAGCAGAG 

R: AGGGCAAGTCATTCCAGTC 

HEX 

59°C 

Diggs and Ardren 
(2008) 

T. arcticus 


Thy62 

(AC)i 7 

F: GGA CGG AGC CAG CAT CAC 

R: GTTT CTT G ATTT C ATA AT C AGGTC A ATAGTC AT 

NED 

59°C 

Susnik et al. (2000) 

T. thymallus 


Tar 106 

(CATA) 24 

F: CGTCCAGTCTGACACAAAG 

R: GTTTCTTATGAACCGAAGGAATCATG 

FAM 

59°C 

Diggs and Ardren 
(2008) 

T. arcticus 


Tth445 

(GATA) 20 

F: TGA CGG CTA CAG GAA TTGT 

R: GTTTCTTCCACAGAGGGTTCTACATTG 

NED 

59°C 

Junge et al. (2010) 

T. thymallus 


Tar 104 

(CATA) 23 

F: TCTTCTCAGTGGCATGACATC 

R: CCTCGTACTCTCTCTTGTCCC 

NED 

60°C 

Diggs and Ardren 
(2008) 

T. arcticus 

5 -Plex 

Thy54 

(TG)i 3 

F: AGAGGGGTCCAGCAACATCA 

R: GTTTGGGGAACCAGTCTAAAGCCT 

HEX 

60°C 

Susnik et al. (1999b) 

T. thymallus 

Tar213 

(CTTT) 33 

F: TTTCCACAGAGGGTTCTACAT 

R: GTTTCTTACTAGAGCAGGGCAGCAGA 

HEX 

60°C 

Diggs and Ardren 
(2008) 

T. arcticus 


Tar414 

(AGAT) 21 

F: GTCGGG AC AT GG ACTCTAC A 

R: GTTTCTTGCAATGCCTCTTTATAGCTT 

NED 

60°C 

Diggs and Ardren 
(2008) 

T. arcticus 


Tth446 

(ATCT) 15 

F: GCCATTCACCCATACTATGC 

R: GTTTCTTCCATTCAGCCACTAGAGC 

FAM 

60°C 

Junge et al. (2010) 

T. thymallus 
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Figure 3. - Median Joining network of haplotypes found in this 
study augmented by additional published haplotypes provided for 
reference. All new haplotypes use the simple abbreviation “Ht” 
with a serial number reflecting the order in which they were first 
found in our data. The small black dots indicate substitutional steps 
whereas the small red dots represent multifurcating nodes. The fol¬ 
lowing haplotypes and corresponding clade names stem from Weiss 
et al. (2002). “ Danube drainage (Northern Alps)”: Dal, Da2, Da4, 
Dali; “Mixed Central Europe ”: Rhl, Rh4, Rh6 all first reported 
from the Rhone basin; “Mixed Central Europe ”: Atl4, At 15 both 
frequent and widespread in the Rhine basin; “Scandinavia”: At6 
reported from Finland; and ATI first reported from the Loire River. 

torial Correspondence Analysis (FCA) with the GENETIX 
software package. Differences among samples are based on 
the presence and absence of alleles, and are displayed in two 
two-dimensional scatter plots using the first four vectors. To 
determine the number of genetically homogenous groups of 
individuals (i.e. clusters) and estimate potential introgres- 
sion of hatchery strains into wild populations we used the 
Bayesian assignment procedure implemented in the program 
Structure v.2.3.3 (Pritchard et al., 2000). The procedure 
assigns individuals according to their genotypes at multiple 
loci into K clusters without using any a priori population 
information. We ran Structure using the admixture-model 
with correlated allele frequencies (Falush et al., 2003). A 
Markov chain Monte Carlo (MCMC) sampling procedure 
with 200,000 iterations after a burn-in of 50,000 iterations 
was used for estimating posterior probabilities. We carried 
out ten replicate runs for each estimation of K ranging from 
1 to 14. We used the method of Evanno et al. (2005), imple¬ 
mented in the program Structure Harvester v.0.6.6 (Earl and 
von Holdt, 2012), to determine the most likely value of K 
(http://taylor0 .biology.ucla.edu/struct_harvest/). 
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Figure 4. - Factorial Correspondence Analyses (FCA) of individu¬ 
als based on the presence and absence of micro satellite alleles. A: 
Bi-variate plot of the first two FCA factors (FI & F2), graphically 
partitioned by population. B: Bi-variate plot of the third and fourth 
FCA factors (F3 & F4) partitioned by population. The percentage 
inertia for each factor is as follows: FC1 = 4.44%; FC2 = 4.03%; 
FC3 = 3.46%; FC4 = 2.82%. 


mtDNA analysis 

The complete mtDNA control region (CR) was amplified 
in 144 individuals (including 34 used in previous reports) 
using the primers LRBT-25 (Forward, 5’-AG A GCG CCG 
GTG TTG TAATC-3’) and LRBT-1195 (reverse, 5’-GCT 
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Figure 5 - Results of the Structure anal¬ 
ysis (Q-values) for 8 loci (A) and 10 
loci (B). Colours represent the different 
“k” units chosen using the approach of 
Evanno et al. (2005). Population names 
are shown on the x-axis, and Q-values 
(contribution or assignment from each 
k partition) are shown on the Y-axis. 
For 8 loci, k = 8, and for 10 loci k= 11. 


AGC GGG ACT TTC TAG GGT C-3’; Uiblein et al ., 2001). 
Amplification was performed with a total volume of 10 pi 
per sample, each consisting of 1 pi diluted DNA extract, 2 pi 
5x Phusion® HF Reaction Buffer, 0.2 pi dNTPs (10 pM), 
0.5 pi of each primer (10 pM.), 0.25 pi Phusion® High 
Fidelity Polymerase (2 units///l) and 5.55 pi of deionized 
water. The amplification temperature profile was as follows: 
after an initial denaturation at 98°C for 30 sec, 35 cycles 
of 10 sec at 98°C, 30 sec at 55°C and 30 sec at 72°C were 
conducted followed by a final elongation step at 72°C for 
10 min and a final hold at 8°C. The sequencing reaction was 
carried out with the internal primers D-Loop-Int5’ (5’-ATA 
TAA GAG AAC GCC CGG CT-3’) for the first part of the 
5’ end, CR-Int-3F-thy (5’- GAA ACC ACT CAC TGAAAG 
CCG-3’) for middle of the PCR product, and CRII-Int-2F- 
thy (5’-GGA ATC CCC CGG CTT CTA C-3’) for the 3’ end. 


After sequencing, DNA fragments were purified with Sepha- 
dexTM G-50 (Amersham Biosciences) following the manu¬ 
facturer’s instructions and visualized on a 3130 XL capillary 
sequencer (Applied Biosystems). Sequences were aligned 
visually using MEGA software v.5 (Tamura et al., 2007). 
Subsequent to alignment, the poly-T region (see Susnik et 
al., 2006) was removed, leaving a total of 1048 bp for anal¬ 
ysis. A median joining haplotype network was constructed 
using the program Pop Art (http://popart.otago.ac.nz). 

RESULTS 

Two of the twelve micro satellites (Tth445 and Tth213) 
were removed from the analysis due to lack of polymor¬ 
phism and problems with null alleles. The number of indi- 
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Table IV. - Occurrence of the different haplotypes in the populations of the Upper Vienne basin, in three other populations of the Loire 
drainage basin, in the Kohlenbach (Danube basin), in the Ain River (Rhone basin), and in four fish farms (in Italic). 



Haplotypes 


Loire basin 

Danube basin 

Rhone basin 

Rhine 



Vienne 

? 

Allier 

Loire 

basin 

Total 

Population 

Ht24 

Ht25 

Ht26 

Ht29 

Ht30 

Ht23 

Ht20 

Ht21 

Atl 

Dal 

Da4 

Dali 

Da30 

Rhl 

Rh4 

Ht4 

Htl4 

Ht3 

Ht27 


Combade 

9 

2 


















11 

Vienne upstream 

11 



















11 

Vienne downstream 

9 


4 

















13 

Maulde 

2 



















2 

Allier Vabres 



3 

4 

3 















10 

Upper Loire 









19 











19 

Ance du N. 







6 

1 












7 

Ain aval 














12 

1 

3 

1 



17 

Besse 






4 






7 








11 

Obenheim 


















3 

8 

11 

Chauvey 










5 

1 


1 







7 

Chazey-Bons 














1 


9 




10 

Kohlenbach 










12 

1 

2 








15 

Total 

31 

2 

7 

4 

3 

4 

6 

1 

19 

17 

2 

9 

1 

13 

1 

12 

1 

3 

8 



viduals successfully processed for the remaining 10 micro¬ 
satellites or mtDNA is listed in table II. For mtDNA, 37 of 
the 87 individuals (43%) sampled from the Upper Vienne 
basin were analysed, which should provide a good estimate 
of haplotype distribution among these populations. Of these 
37 individuals, only three different haplotypes, one already 
observed in the Allier River in 2011 (Ht26, Persat et al., 
2013), and two unknown up to now (Ht24 and Ht25) were 
revealed (Tab. IV). Among the latter, one is common in all 
the Vienne populations whereas the second appeared only at 
low frequencies in the Combade River. None of these hap¬ 
lotypes appeared among the nine haplotypes recorded from 
the four fish farms. While eight of the nine fish-farm hap¬ 
lotypes grouped with clades not representative of the Loire 
(i.e. Danube, Rhone and Rhine basins), one newly reported 
Loire-clade haplotype (Ht23) was found in four individuals 
of the Besse fish farm. 

Based on the median-joining nework (Fig. 3), the Loire 
haplotypes form a very distinct clade, and this clade is well 
structured reflecting each of the three main Loire sub-drain¬ 
ages still inhabited by grayling. Thus, the two new “Vienne” 
haplotypes (Ht24, Ht25) form a distinct branch within the 
Loire clade, whereas the third “Vienne” haplotype (Ht26) 
is closely related to the other Allier branch haplotypes and 
the Loire clade haplotype from the Besse fish farm (Ht23) 
as well. In summary, at the mtDNA level, there was no evi¬ 
dence of allochthonous introgression into the Upper Vienne 
basin. 

Based on microsatellite data, a number of samples devi¬ 
ated from HWE with positive Fis-values, including the 
Combade sample, and two of the four fish farms (Besse and 


Chazey-Bons) suggesting substructure, whereby the Besse 
fish farm sample additionally had one locus with a signifi¬ 
cant negative Fjs-value (Tab. V). In contrast, the upper Loire 
sample significantly deviated from HWE with a negative 
F\ s-value suggesting, together with the low allelic diver¬ 
sity, family structure within the sample. The mean number 
of alleles for Loire basin samples (i.e. Upper Vienne, Upper 
Loire and Allier) was rather low compared to the other sam¬ 
ples, but within the Loire basin, the Vienne reach appeared 
to be more diverse than the Upper Loire (Tab. V). The high¬ 
est expected heterozygosity was found in the Besse fish-farm 
sample, being higher than the most diverse natural popula¬ 
tion in this study, the Lower Ain River. However, consider¬ 
ing the deviation from HWE, this sample is likely composed 
of more than one source population. 

The lack of variation at the Thyl seems to be specific to 
the Upper Vienne drainage basin, while the low variability 
of Thy54 and Thtl23, corresponds to their usual low vari¬ 
ability elsewhere. 

The Factorial Correspondence Analysis based on the 
presence and absence of alleles within individuals provides 
an overview of the genetic relationships among all sam¬ 
ples (Fig. 4A, B). The first factor strongly discriminates 
individuals from the Loire basin (Vienne, Upper Loire and 
Upper Allier) from all others. The second factor especially 
discriminates Danubian samples from Rhine and Rhone 
(i.e. Ain, Chazey-Bons, Obenheim), whereas the third fac¬ 
tor distinguishes the Ain and Chazey-Bons (both of Rhone 
origin) from the Obenheim (Rhine origin) samples. Finally, 
the fourth factor discriminates the Loire samples from each 
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Table V. - Summary statistics of the 10 micro satellite loci for each population. N: individual number of samples; A N : number of alleles; 
A r : allelic range; H 0 \ observed heterozygosity; H E : expected heterozygosity; F IS : inbreeding. The Upper Vienne sample sites are listed both 
separately and as one pooled group. Bold numbers with an asterisk indicate a statistically significant deviation from HWE. 



N 


TarlOl 

Tar 106 

Thyl 

Thy62 

Tth305 

Tth447 

Tar 104 

Thy54 

Tth414 

Tth446 

Overall 

Upper 

37 

An 

1 

8 

1 

5 

6 

4 

9 

3 

1 

5 

55 

Vienne 


Ar 

260-336 

191-307 

109-109 

202-256 

274-316 

169-179 

173-241 

191-199 

356-502 

227-251 


(Combade) 


Ho 

0.737 

0.583 

0 

0.568 

0.546 

0.000 

0.733 

0.629 

0.565 

0.639 

0.554 



He 

0.748 

0.682 

0 

0.562 

0.632 

0.540 

0.668 

0.613 

0.686 

0.745 

0.217 



Fis 

0.042 

0.159 

NA 

0.003 

0.160 

0.132 

-0.080 

-0.011 

0.198 

0.156 

0.087* 

Upper 

34 

A n 

3 

6 

1 

4 

6 

3 

6 

3 

10 

5 

47 

Vienne 


Ar 

260-284 

191-307 

109-109 

204-256 

274-316 

169-179 

185-241 

191-199 

356-514 

227-251 


(Vienne 


H 0 

0.450 

0.677 

0 

0.677 

0.552 

0.618 

0.765 

0.647 

0.680 

0.765 

0.583 

upstream) 


H e 

0.405 

0.751 

0 

0.653 

0.617 

0.629 

0.731 

0.635 

0.770 

0.695 

0.589 



Fis 

-0.086 

0.114 

NA 

-0.022 

0.123 

0.033 

-0.032 

-0.003 

0.137 

-0.082 

0.026 

Upper 

15 

An 

5 

6 

1 

4 

4 

3 

6 

3 

7 

5 

44 

Vienne 


Ar 

260-336 

191-307 

109-109 

204-256 

274-316 

169-179 

173-241 

191-199 

356-502 

227-251 


(Vienne 


H 0 

0.615 

0.867 

0 

0.800 

0.692 

0.800 

0.583 

0.357 

0.667 

0.667 

0.605 

downstream) 


He 

0.701 

0.738 

0 

0.651 

0.660 

0.640 

0.618 

0.554 

0.821 

0.738 

0.612 



Fis 

0.162 

-0.141 

NA 

-0.196 

-0.009 

-0.217 

0.099 

0.387 

0.244 

0.130 

0.053 

Upper 

2 

An 

1 

2 

1 

1 

1 

1 

2 

3 

1 

2 

15 

Vienne 


Ar 

276-276 

191-231 

109-109 

256-256 

304-304 

175-175 

185-225 

191-199 

502-502 

239-251 


(Maulde) 


H 0 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 



H e 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 



Fis 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

Upper 

88 

An 

7 

9 

1 

6 

1 

4 

11 

4 

13 

5 

67 

Vienne 


Ar 

260-336 

191-311 

109 

202-256 

274-316 

169-179 

173-241 

191-199 

356-514 

227-276 


(sub-total) 


Ho 

0.574 

0.667 

NA 

0.636 

0.561 

0.602 

0.718 

0.600 

0.610 

0.690 

0.629 



H e 

0.672 

0.7544 

NA 

0.647 

0.632 

0.649 

0.737 

0.644 

0.816 

0.767 

0.702 



Fis 

0.154* 

0.122* 

NA 

0.023 

0.120 

0.161 

-0.033 

0.074 

0.260* 

0.106* 

0.111* 

Upper Loire 

43 

An 

5 

8 

3 

4 

5 

4 

4 

3 

5 

4 

45 



a 

Grange 

256-348 

191-307 

109-125 

204-264 

274-316 

169-185 

185-225 

191-197 

350-478 

235-243 




Ho 

0.615 

0.698 

0.372 

0.762 

0.7692 

0.767 

0.707 

0.721 

0.442 

0.721 

0.657 



H e 

0.649 

0.643 

0.391 

0.717 

0.726 

0.617 

0.583 

0.546 

0.371 

0.623 

0.586 



Fis 

0.065 

-0.074 

0.059 

-0.051 

-0.046 

-0.232* 

-0.202* 

-0.310* 

-0.180 

-0.146 

-0.109* 

Allier 

42 

An 

7 

7 

3 

5 

5 

5 

6 

4 

11 

4 

57 



Ar 

252-296 

191-311 

109-129 

230-260 

284-316 

169-181 

169-241 

189-199 

356-514 

227-243 




Ho 

0.600 

0.825 

0.683 

0.725 

0.775 

0.590 

0.564 

0.100 

0.821 

0.575 

0.626 



H e 

0.618 

0.793 

0.558 

0.710 

0.738 

0.530 

0.599 

0.120 

0.818 

0.575 

0.606 



Fis 

0.042 

-0.027 

-0.213* 

-0.009 

-0.037 

-0.033 

0.071 

0.177 

-0.010 

0.013 

-0.020 

Ain 

43 

An 

16 

10 

10 

4 

14 

7 

19 

11 

18 

8 

117 



Ar 

208-324 

191-299 

113-157 

188-208 

282-350 

160-190 

157-285 

189-211 

366-558 

215-259 




H 0 

0.850 

0.714 

0.781 

0.325 

0.825 

0.275 

0.805 

0.905 

0.691 

0.738 

0.691 



H e 

0.882 

0.765 

0.795 

0.326 

0.777 

0.251 

0.891 

0.845 

0.757 

0.801 

0.709 



Fis 

0.049 

0.079 

0.031 

0.016 

-0.049 

-0.08 

0.109* 

-0.058 

0.100 

0.091 

0.038 

Besse 

29 

An 

7 

12 

7 

9 

7 

8 

10 

13 

9 

6 

88 



Ar 

204-324 

239-283 

109-136 

202-234 

284-382 

156-176 

149-231 

177-231 

348-494 

219-279 




Ho 

0.593 

0.966 

0.828 

0.857 

0.546 

0.655 

0.552 

1.000 

0.417 

0.370 

0.678 



H e 

0.723 

0.858 

0.823 

0.760 

0.778 

0.725 

0.876 

0.885 

0.789 

0.630 

0.785 



Fis 

0.198* 

-0.108 

0.012 

-0.110 

0.320* 

0.038 

0.386* 

-0.112* 

0.488* 

0.428* 

0.154* 
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Table V. - Continued. 



N 


TarlOl 

Tar 106 

Thyl 

Thy62 

Tth305 

Tth447 

Tar 104 

Thy54 

Tth414 

Tth446 

Overall 

Chauvey 

33 

A n 

12 

6 

5 

7 

0 

3 

1 

11 

3 

0 

54 



Ar 

264-352 

175-283 

110-134 

214-248 

NA 

152-178 

149-239 

175-235 

348-356 

NA 




H 0 

0.833 

0.781 

0.697 

0.656 

NA 

0.424 

0.394 

0.844 

0.505 

NA 

0.642 



H e 

0.840 

0.709 

0.702 

0.585 

NA 

0.574 

0.346 

0.828 

0.548 

NA 

0.642 



Fis 

0.025 

-0.086 

0.022 

-0.105 

NA 

0.038* 

-0.124* 

-0.004 

-0.069 

NA 

-0.002 

Obenheim 

40 

A n 

3 

10 

4 

5 

8 

7 

12 

10 

12 

9 

80 



Ar 

208-320 

187-215 

113-145 

202-214 

366-414 

152-196 

149-285 

181-235 

352-526 

207-275 




Ho 

0.275 

0.950 

0.275 

0.575 

0.889 

0.825 

0.923 

0.790 

0.526 

0.872 

0.690 



H e 

0.247 

0.849 

0.248 

0.660 

0.804 

0.759 

0.802 

0.708 

0.813 

0.843 

0.673 



Fis 

-0.103 

-0.106* 

-0.097 

0.141 

-0.091 

-0.073 

-0.139* 

-0.102 

0.364* 

-0.022 

-0.012 

Chazey-Bons 

43 

An 

10 

13 

10 

8 

0 

8 

14 

10 

13 

6 

92 



Ar 

196-316 

215-291 

110-157 

188-230 

NA 

156-188 

149-273 

177-231 

348-538 

227-279 




Ho 

0.564 

0.875 

0.585 

0.595 

NA 

0.644 

0.884 

0.814 

0.707 

0.368 

0.671 



H e 

0.823 

0.864 

0.806 

0.620 

NA 

0.605 

0.868 

0.770 

0.829 

0.741 

0.769 



Fis 

0326* 

0 

0.285* 

0.051 

NA 

-0.073 

-0.007 

-0.046 

0.159* 

0.512* 

0.151* 

Kohlenbach 

25 

A n 

14 

8 

5 

6 

0 

4 

8 

8 

4 

0 

57 



Ar 

208-352 

195-291 

110-140 

202-224 

NA 

152-170 

149-255 

175-229 

348-386 

NA 




Ho 

0.960 

0.840 

0.708 

0.625 

NA 

0.667 

0.739 

0.826 

0.636 

NA 

0.750 



H e 

0.894 

0.825 

0.626 

0.705 

NA 

0.561 

0.682 

0.771 

0.656 

NA 

0.715 



Fis 

-0.054 

0.002 

- 0.111 

0.134 

NA 

-0.207 

-0.106 

-0.049 

0.053 

NA 

-0.028 


Table VI. - Pairwise Fst distances in eight populations. Significant F$t values are indicated in bold letters. 



Upper Loire 

Allier 

Ain 

Besse 

Chauvey 

Obenheim 

Chazey-Bons 

Kohlenbach 

Upper Vienne 

0.261 

0.260 

0.285 

0.243 

0.356 

0.339 

0.265 

0.320 

Upper Loire 


0.312 

0323 

0.265 

0356 

0361 

0.302 

0.339 

Allier 



0.308 

0.225 

0373 

0341 

0.283 

0.334 

Ain 




0.204 

0324 

0.264 

0.085 

0.284 

Besse 





0.210 

0.179 

0.136 

0.172 

Chauvey 






0342 

0.252 

0.119 

Obenheim 







0.240 

0.296 

Chazey-Bons 








0.209 


other, primarily Allier from Upper Loire with the Upper 
Vienne locations in an intermediate position. 

Two series of Structure analyses were carried out, one 
with 10 loci, and one with the loci Thy 1 and Tth305 removed 
due to their lack of polymorphism in several populations. 
For the reduced data set, the Evanno approach revealed 
K = 8 as the most likely number of populations (Fig. 5A). 
Hereby the various samples sites of the Upper Vienne sys¬ 
tem appear as one homogeneous group, and the other two 
sub-drainages of the Loire catchment (Upper Loire and 
Allier) each appear separately as homogeneous groups. 
Collectively, these results strongly support the absence of 
non-Loire genotypes throughout the Loire catchment. Two 
additional genetic groups are clearly revealed representing 
the Obenheim hatchery, and the Danubian Kohlenbach pop¬ 
ulation together with the Chauvey-bons hatchery confirming 


the notion that this hatchery stock was founded with Danube 
drainage grayling. Additional genetic groups were identified 
in the samples from Chazey-Bons and Besse hatcheries but 
these populations are clearly made up of multiple genotype 
groups. Finally, the Ain population, while revealing an addi¬ 
tional genetic subdivision in this Structure analysis, shows 
admixture from genotypes apparently common to those in 
the Chazey-Bons hatchery. Re-running the structure with all 
10 loci merely shows additional resolution with the genetic 
subdivisions of the Besse and Chazey-Bons hatcheries with 
no effect on the genetic composition or inferences concern¬ 
ing the Loire catchment (Fig. 5B). Overall, moderate to high 
Fst -values for nearly all pairwise comparisons (Tab. VI) also 
supported the genetic differences seen among populations in 
the Structure analysis. 
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DISCUSSION 

The genetic analysis of four European grayling popula¬ 
tions from the Upper Vienne River drainage basin using both 
nDNA (microsatellites) and mtDNA markers shows that they 
all share a typical Loire (ligerian) genetic profile, but with 
sufficient differentiation from other sub-drainages to strong¬ 
ly support that they represent pure or nearly pure native 
stocks. Despite decades of stocking, much of which likely 
implied foreign strains of grayling, no sign of introgres- 
sion from foreign sources was detected. The occurrence of 
a shared haplotype between the lowermost population in the 
Vienne River and fish from the Allier River is likely based 
on human-induced exchanges within the Loire catchment 
(e.g. from the Besse farm), but could also be explained by 
historical exchange or retained ancestral polymorphisms. We 
note that these river stretches are separated by at least 1000 
km of river, most of which is heavily dammed so that gene 
exchange in modern times is virtually impossible. While the 
above-noted haplotype was not found in our sample from 
the Besse population, the sample from this hatchery need 
not be representative of the whole production of the facil¬ 
ity. The Besse hatchery lineage has a long and complex his¬ 
tory stemming back to the early rearing attempts at the CSP 
fish farm of Augerolles (Carmie et al., 1985; Persat et al ., 
1997), which mainly started with imported eggs before shift¬ 
ing to local spawners from the Loire basin. Our current data 
set, however, is likely insufficient to evaluate quantitatively 
the lineage composition of the Besse hatchery. Nonetheless, 
it is clear that the current sample from the Besse hatchery 
reveals at least two genotypic groups that do not appear in 
any sample from the Loire catchment. Likewise, none of the 
other three hatcheries, each of which is genetically distinct, 
appears to have contributed to the current genetic architec¬ 
ture of the upper Loire. 

The distinct genotypic group in the Structure analysis 
combined with the mtDNA haplotype data strongly support 
that the Obenheim hatchery consists of Rhine basin grayling, 
and this confirms the report that the strain originated from 
wild-caught spawners from the Steina brook, a tributary of 
the Wutach River in Baden-Wiirttemberg (Martin Gerber, 
pers. com.). Thus, we have highly recognizable genotypes 
(microsatellites) and haplotypes (mtDNA) from three dis¬ 
tinct basins (Rhone, Rhine and Danube) from outside of the 
Loire catchment. 

There is less than 100 km of grayling habitat currently 
in the Upper Vienna region (Lig. 1) and thus stocking den¬ 
sities roughly reached up to > 1000 individuals per km in 
some years (Tab. I), for the entire region. Densities of wild 
fish are not available, but are considered low. Overall, the 
most likely explanation for the lack of a genetic trace for the 
largely non-native strains of stocked fish is very low survival 
after release. 


The mechanisms responsible for the low survival are not 
evaluated in this study, but more than likely reflect the typical 
suite of factors revealed for other salmonid species. That is, 
a combination of maladaptation (genetic) for the given envi¬ 
ronment, behavioural or physiological imprinting in captiv¬ 
ity (e.g. adaption to artificial diet, still versus running water), 
and poor competitiveness with the existing native stocks. 
Some observations on the Rhone River by-passed section in 
Chautagne (Savoy; Persat et al ., 1997) shown that of 54 000 
fingerlings stocked in late August 1993, an unknown propor¬ 
tion was still present three months later despite of a severe 
period of high water levels soon after having being released. 
However, they were much smaller sized (9.8 cm versus 
14.9 cm) and occupied a much more quiet habitat (shal¬ 
lows at 10 cm depth with a rather low water velocity) along 
the gravel shore compared to typical habitat of the resident 
native grayling that are found in the main riffles of the river 
at 40-80 cm depth (Mallet et al., 2000). Lrom our study on 
the Rhone in Chautagne, we assume that the stocked gray¬ 
ling occupied sub-optimal to poor habitat and thus were 
marginalized and experienced very poor survival (no trace 
of survival of the fishes was found in the years thereafter). 

While from a conservation standpoint it is positive that 
genetic introgression has not taken place, we should recog¬ 
nize that this does not exclude the possibility that continuous 
stocking has had negative consequences in raising short-term 
intraspecific competition, introducing pathogens, promoting 
higher densities of predators, and elevating angler expecta¬ 
tions. 

Ongoing investigations on other grayling populations 
in Lrance seem to confirm that the competitive advantage 
of native stocks over translocated ones is often the rule; 
observed examples of successful stocking programs are 
mainly found in rivers free of native populations such as 
in the Seine and Garonne basins (Persat, 2011). However, 
wide-scale contact and interbreeding is reported between 
close lineages in the upper Rhone and Ain river basins where 
Rhone related populations more or less mix with genotypes 
of Swiss-Rhine origin, especially those found in the Geneva 
Lake drainage basin (Cattaneo et al., 2011). Both Cattaneo 
et al. (2011) and Vonlanthen and Salzburger (2011), with a 
more detailed sampling, also clearly demonstrated the sur¬ 
vival of stocked Rhone basin grayling in the Birs River of the 
Rhine catchment in Switzerland, though the river is severe¬ 
ly fragmented by dams, and not all stretches demonstrated 
the sympatry of native and stocked lineages. In contrast, 
both Cattaneo et al. (2011) and Vonlanthen and Schlunke 
(2015) reported no genetic trace of multiple stockings from 
multiple sources (for the latter study, three hatcheries were 
documented) throughout the Aar, Limmat and Reuss rivers 
in the Rhine catchment of Switzerland. Lor the latter study, 
the documented stocked fish were of the same major line¬ 
age (Rhine catchment), but nonetheless clearly distinguish- 
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able with microsatellite markers. Stocking densities ranged 
from 30 000 to over 120 000 per year for 18 years, across 
river reaches of ca. 100 km. Yet, as it was the case in our 
study, no trace of the hatchery material was found. Thus 
despite knowledge that stocked grayling of various diver¬ 
gent origins do introgress or even in some situations replace 
native stocks (e.g. Meraner et al., 2014; Weiss et al., 2015) 
our study and the cited reports provide multiple examples 
of complete failure of stocked lineages to survive and intro¬ 
gress into native populations. Currently, we lack insight into 
the decisive mechanisms that support or hinder survival and 
introgression of stocked grayling. Nonetheless, a frequently 
cited factor for the promotion of hybridization and introgres¬ 
sion, and one that also has at least some supportive data for 
salmonid fishes is the environmental conditions and demo¬ 
graphic health of the native populations (Rhymer and Sim- 
berloff, 1996; Seehausen et al., 2008; Bittner et al., 2010; 
Winkler et al., 2011). That is, the less disturbed systems with 
healthy native populations present a more likely barrier to 
introgression than disturbed environments with low popula¬ 
tion density. 

The efficacy of management focused on stocking hatch¬ 
ery grayling is thus called into question, especially for flow¬ 
ing water habitats that contain native stocks. This once again 
calls for alternative investment in improved habitat condi¬ 
tions, including restoration of the riverbed and connectiv¬ 
ity, reduction of hydropeaking, increase of residual flows 
and improved monitoring and control of reservoir flushing 
operations. 

Finally, if a stocking appears definitely necessary in a 
reach where the grayling is obviously extinct, it should be 
based on translocation of wild fry as young as possible (such 
as at emergence) because Bardonnet and Gaudin (1990) 
showed that even the artificial conditions during the incuba¬ 
tion could alter the behaviour and emergence timing of the 
larvae. Of course, these fry should be collected in the same 
catchment, and as close as possible to the reach to stock in 
order to have the highest likelihood of introducing geno¬ 
types that are adapted to the local environmental conditions. 
Another promising approach that is gaining popularity in 
central Europe is referred to as “cocooning”. The approach is 
based on an adaptation of so-called Vibert or Whitlock-Vib- 
ert breeding boxes (Vibert, 1953; Whitlock, 1978), where¬ 
by several changes to these designs as well as the general 
application strategy have proved to be very successful for 
several salmonid species including grayling (Holzer et al., 
2011). The overall concept supports the maximal degree of 
adaptation and conditioning to the local environment as the 
fish are hatched in the riverbed itself with the only artificial 
component of the procedure being the choice and stripping 
of parental fish. 
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